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LONG-TERM  GOAL 

Synthetic  Aperture  Radar  (SAR)  imagery  has  shown  great  promise  for  depicting  the  vast  array  of 
phenomena  that  govern  the  behavior  of  the  ocean  mixed  layer  and  marine  atmospheric  boundary 
layer  (MABL).  (Alpers,  et  al.,  1981;  Beal,  et  al.,  1981;  Vesecky  and  Stewart,  1982).  The 
variation  of  the  backscattered  intensity  field  depicted  in  SAR  imagery  is  directly  related  to  the 
horizontal  distribution  of  those  sea-surface  roughness  elements  having  scales  generally 
comparable  to  the  wavelength  of  the  radiation  transmitted  by  the  SAR.  The  local  amplitude  of 
the  (centimeter-scale  for  most  SAR  systems)  surface  waves  that  produce  this  roughness  depends 
on  a  broad  range  of  oceanic  and  atmospheric  processes  and  their  interactions  (e.g.,  Elachi,  1987). 
Because  these  waves  are  driven  by  the  surface  stress  and  locally  modulated  by  wave-current 
interactions  and  surfactant  slicks,  SAR  images  frequently  reveal  features  related  to 
oceanographic  processes  such  as  current  boundaries,  internal  waves,  or  tidal  flow  over 
bathymetry,  as  well  as  variations  in  the  surface  stress  due  to  atmospheric  processes. 

Our  long-term  goal  in  this  research  effort  is  to  utilize  the  multiscale  information  in  the 
atmospheric  signatures  on  SAR  images  to  diagnose  a  quantitative  description  of  the  MABL, 
including  the  depth,  stability,  wind  speed,  wind  direction,  sea-surface  stress,  and  buoyancy  flux 
on  the  mesoscale  and  submesoscale.  Because  of  its  potential  for  yielding  both  boundary  layer 
depth  and  the  surface  wind  field  at  high  horizontal  resolution,  this  application  of  SAR  data 
represents  a  significant  and  innovative  advance  over  most  scatterometer  algorithms  that  yield 
only  coarse-resolution  wind  fields.  Moreover,  because  conventional  scatterometry  cannot  resolve 
the  turbulence  structures  in  the  MABL,  it  cannot  be  used  to  diagnose  the  surface  layer  stability 
and  so  cannot  yield  wind  speed  estimates  corrected  for  this  important  effect. 

SCIENTIFIC  OBJECTIVES 

SAR  has  the  potential  to  overcome  some  of  the  inherent  limitations  of  conventional 
scatterometry  by  providing  information  on  the  kilometer-scale  variability  of  surface  stress.  This 
variability  is  directly  related  to  the  intensity  of  the  primary  turbulence  structures  in  the 
convective  MABL  via  well-known  similarity  relationships  (e.g.,  Panofsky  and  Dutton,  1984; 
Stull,  1988).  These  structures,  exemplified  by  three-dimensional  convective  cells  and  two- 
dimensional  longitudinal  rolls  (Woodcock,  1975)  strongly  modulate  the  sea-surface  stress  and  so 
leave  SAR-detectable  footprints  on  the  sea  surface  in  the  form  of  kilometer-scale  patterns.  By 
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quantitatively  analyzing  the  variability  in  SAR  backscatter  intensity  wrought  by  these  structures, 
we  are  developing  methods  for  determination  of  boundary  layer  depth,  boundary  layer  stability, 
boundary  layer  wind  direction,  and  surface  layer  wind  speed. 

APPROACH 

We  obtain  the  boundary  layer  depth  directly  from  the  dominant  horizontal  scale  of  the  SAR 
backscatter  variability  using  a  similarity  relationship  for  the  aspect  ratio  of  the  horizontal 
wavelength  to  boundary  layer  depth  of  the  associated  three-dimensional  convective  cells  (Sikora, 
1996a;  Sikora,  1996b;  Sikora  and  Young,  1996;  Sikora,  et  al.,  1997).  We  obtain  the  wind 
direction  by  noting  that  the  footprints  of  the  MABL  coherent  structures  are  symmetric  in  the 
direction  perpendicular  to  the  mean  wind,  but  asymmetric  along  the  mean  wind  (Sikora,  et  al. 
1995).  We  use  two-dimensional  (2D)  spectral  analysis  to  extract  both  the  dominant  cell 
horizontal  wavelength  and  the  dominant  direction  of  the  features. 

A  primary  focus  of  our  research  is  the  development  of  an  iterative  algorithm  using  Monin- 
Obukhov  similarity  relationships  to  relate  surface  stress  estimates  to  surface  layer  wind  speed. 
Monin-Obukhov  similarity  captures  the  effect  of  surface  layer  stability  on  the  ratio  of  boundary 
layer  turbulence  intensity  to  mean  wind  speed.  Thus,  it  is  an  appropriate  tool  for  relating  the 
degree  of  kilometer-scale  variability  in  the  SAR  backscatter  intensity  to  the  mean  backscatter 
intensity.  By  inverting  this  relationship,  we  can  diagnose  the  surface  layer  stability  (Monin- 
Obukhov  length)  from  the  ratio  of  the  first  and  second  moments  of  the  SAR  backscatter 
distribution.  An  iterative  solution  of  similar  systems  of  Monin-Obukhov  similarity  relations  has 
been  done  for  a  number  of  applications  (e.g.,  Young  and  Kristensen,  1992;  Fairall,  et  al.,  1996). 
As  with  other  algorithms  of  this  ilk,  ours  starts  with  a  first  guess  that  is  a  semi-reasonable  value 
for  the  surface  layer  stability.  Convergence  of  the  method  is  not  dependent  upon  our  having  a 
particularly  accurate  first  guess  of  this  parameter;  both  neutral  stability  and  the  climatological 
average  stability  value  have  been  used  successfully  in  previous  studies. 

We  are  using  the  nonlinear  convective  structure  model  (Shirer,  et  al.,  1995,  1996,  1997)  primarily 
to  provide  the  theoretical  underpinnings  for  the  refinements  of  current  atmospheric  boundary 
layer  similarity  theory  needed  to  refine  our  new  iterative  algorithm.  This  task  requires  that  we 
explore  as  thoroughly  as  possible  the  relevant  parameter  space  (given  for  example  by  air- sea 
temperature  difference  and  surface  wind  speed).  Our  highly  truncated  spectral  model,  developed 
under  previous  ONR  funding  (Shirer,  et  al.,  1995,  1996),  resolves  only  the  dominant  (kilometer- 
scale)  boundary  layer  eddies,  and  so  fills  this  role.  The  essential  surface  layer  physics  is 
incorporated  using  Monin-Obukhov  similarity-based  lower  boundary  conditions  linking  the 
model-resolved  flow  to  surface  fluxes  of  momentum  and  heat  (Zuccarello,  1994).  As 
demonstrated  by  Lambert  (1995),  the  model  has  the  ability  to  produce  convectively  induced 
kilometer-scale  stress  variability  patterns  resembling  those  seen  on  the  SAR  image  studied  by 
Sikora,  et  al.,  (1995).  We  employ  linear  analysis  of  the  model  solutions  to  obtain  the  preferred 
wavelengths  and  orientations  of  the  axes  of  symmetry  in  the  stress  patterns  for  use  in  improving 
the  similarity  theory  results  as  described  above.  We  then  employ  the  nonlinear  solutions 
themselves  to  refine  the  rstability,  and  air/sea  buoyancy  flux  from  SAR  imagery  of  the  sea 
surface  under  unstable  conditions.  A  detailed  description  of  this  algorithm  as  well  as  results  of  the 
sensitivity  analysis  were  reported  at  the  4th  International  Conference  on  Remote  Sensing  for 
Marine  and  Coastal  Environments  (Young,  et  al.,  1997a).  Three  papers  were  presented  at  the 
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American  Meteorological  Society’s  12th  Symposium  on  Boundary  Layers  and  Turbulence.  The 
first  (Young,  et  al.,  1997b)  discussed  postdoctoral  fellow  Todd  Sikora’s  preliminary  testing  of  this 
algorithm  on  a  small  suite  of  cases  from  HI  RES  II.  The  second  (Shirer,  et  al.,  1997)  summarized 
graduate  student  David  Beberwyk’s  case  study,  using  the  convective  structures  model,  of  an 
ERS-1  overpass  of  the  northern  part  of  the  Gulf  Stream  during  the  HI  RES  II  experiment.  The 
third  (Winstead,  et  al.,  1997)  reviewed  ASSERT  graduate  student  Nathaniel  Winstead’s 
observational  study  of  the  SAR  signatures  of  drainage-flow-induced  exit  jets  over  Chesapeake 
Bay. 

RESULTS 

George  Young  and  Todd  Sikora  developed  a  similarity-based  iterative  solution  technique  for 
deriving  both  the  surface  layer  stability  and  the  stability-corrected  scatterometer  wind  speed  from 
SAR  imagery.  The  sensitivity  of  this  algorithm  to  uncertainty  in  the  SAR-based  estimates  of  the 
boundary  layer  depth  was  evaluated  and  found  to  be  negligible  for  the  scatterometer  winds,  with 
typical  uncertainty  in  the  boundary  layer  depth  leading  to  wind  differences  on  the  order  of  one 
percent.  Thus,  the  algorithm  has  the  capability  of  correcting  scatterometer  winds  for  stability 
effects  even  in  the  face  of  uncertainty  in  the  boundary  layer  depth.  The  surface  buoyancy  flux 
estimate  that  is  also  produced  by  this  algorithm  is  rather  more  sensitive,  displaying  fractional 
uncertainties  only  slightly  less  then  those  of  the  input  boundary  layer  depths.  For  the  magnitude 
of  boundary  layer  depth  uncertainty  observed  in  our  preliminary  tests  at  sea,  this  sensitivity 
implies  air/sea  buoyancy  fluxes  having  levels  of  uncertainty  that  are  comparable  with  those  from 
the  best  bulk  aerodynamic  algorithms.  Unlike  the  bulk  aerodynamic  algorithms,  however,  our 
SAR-based  algorithm  requires  no  in  situ  data  and  so  may  be  used  where  no  surface  platform  or 
low-level  aircraft  can  be  deployed.  Head-to-head  testing  of  our  SAR-based  algorithm  against  in 
situ  wind  and  flux  measurements  began  with  two  ERS-1  images  of  the  northwest  edge  of  the  Gulf 
Stream  from  the  HI  RES  II  experiment.  The  results  from  these  preliminary  test  cases  show  that 
the  stability  correction  to  the  scatterometer-diagnosed  wind  speed  can  be  significant  in  the 
unstable  MABL.  The  algorithm  correctly  detected  the  atmospheric  surface  layer  stability 
transition  at  the  edge  of  the  Gulf  Stream.  Scale  limits  on  the  variance  analyzed  are  needed  to 
eliminate  occasional  sign  errors  in  the  atmospheric  surface  layer  stability  and  should  improve  the 
quantitative  accuracy  of  this  parameter.  The  surface  buoyancy  flux  maps  diagnosed  from  these 
SAR  images  also  captured  the  effect  of  the  Gulf  Stream,  but  their  values  tended  to  be 
systematically  low,  suggesting  again  that  filtering  of  the  SAR  imagery  to  eliminate  non-turbulence 
scales  of  variability  is  required.  These  results  are  expected  because  the  algorithm  uses  mixed- 
layer  similarity  for  diagnosing  atmospheric  characteristics  from  variability  in  the  SAR  imagery. 
Our  next  task  is  to  incorporate  such  filtering  into  the  algorithm  and  to  test  it  against  a  larger  set  of 
images. 

David  Beberwyk  completed  the  convective  structures  model  development  and  performed  a  case 
study  of  the  convective  MABL  observed  during  HI-RES  II.  As  model  input  he  used  data  taken 
from  the  RV  Columbus  Iselin  and  the  first  of  two  radiosondes  launched  from  the  ship.  Spectral 
analysis  of  a  nearly  coincident  ERS-1  SAR  image  by  Todd  Sikora  provided  the  cellular  spacing 
that  the  model  should  give.  The  model  was  able  to  reproduce  this  spacing  for  a  plausible 
boundary  layer  wind  profile,  which  had  to  be  chosen  because  none  was  measured  by  the 
radiosonde,  and  nonlinear  solutions  provided  estimates  of  the  stress  variability  produced  by  the 
MABL  large  eddies  (Beberwyk,  1997;  Shirer,  et  al.,  1997).  The  resulting  stress  patterns  were 
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oriented  within  10  degrees  of  the  10  m  current-relative  wind  in  this  case,  which  matched  the 
observations  very  well. 

AASERT  student  Nathaniel  Winstead  studied  an  ERS-1  SAR  image  taken  late  in  the  evening  of 
May  9,  1992.  This  image  shows  fan-shaped  fingers  of  brightness  flowing  from  the  west  shore  of 
the  Chesapeake  Bay.  These  SAR  signatures  are  perfectly  aligned  with  the  creek  basins  and 
narrow  canyons  present  along  the  shore,  strongly  suggesting  a  link  between  the  topography 
around  the  Bay  and  the  SAR  signatures.  This  coupled  with  the  time  of  the  image  suggests  that  the 
patterns  evident  on  the  SAR  are  the  signatures  of  exit  jets  induced  by  nocturnal  drainage  flow.  A 
regression  analysis  was  performed  to  link  the  size  and  shape  of  the  individual  basins  to  the  length 
of  the  SAR  signatures  (Winstead,  et  al.,  1997).  This  regression  analysis  indicates  that  longer 
basins  with  wider  gap  widths  correspond  to  longer  jets.  Concurrent  with  this  statistical  study,  a 
numerical  model  was  developed  to  study  the  dynamics  of  these  drainage  flow  structures.  The 
ultimate  goal  of  the  research  is  to  develop  a  similarity  theory  linking  the  length  of  the  exit  jets 
with  the  various  basin  characteristics. 

IMPACT/APPLICATION 

Our  new  quantitative  approach  will  not  only  provide  more  accurate  wind  mapping  on  a  higher 
resolution  than  do  existing  scatterometer  algorithms,  but  it  will  also  allow  diagnoses  to  be 
performed  much  closer  to  strong  discontinuities  such  as  coasts  and  ocean  current  boundaries. 
Thus  our  method  offers  the  promise  of  mapping  winds  and  boundary  layer  depth  variations 
within  the  mesoscale  circulations  caused  by  these  surface  discontinuities.  Wind  fields  as  well  as 
maps  of  boundary  layer  depth  and  stability  resulting  from  our  studies  are  expected  to  yield 
interesting  insights  into  the  mesoscale  flow  fields  near  coasts  and  ocean  current  boundaries. 
There  is  a  growing  realization  from  experiments  such  as  HI-RES  (Sublette,  1994)  that  current 
boundaries  lead  to  mesoscale  solenoidal  wind  circulations  and  frontal  systems  just  as  significant 
as  those  associated  with  coasts. 

TRANSITIONS 

We  foresee  a  number  of  operational  uses  of  our  results.  First,  they  will  help  improve  the 
initialization  of  horizontally  varying  fields  of  boundary  layer  depth,  surface  wind,  and  air/sea 
temperature  difference  in  operational  numerical  weather  prediction  models,  especially  in  data- 
sparse  regions.  Second,  the  high-resolution  wind  fields  could  be  used  to  correct  synoptic -sc ale 
fields  that  drive  global  ocean  wave  models.  Third,  these  high-resolution  wind  estimates  would 
provide  important  constraints  on  coastal  remote  sensing  applications  such  as  the  imaging  of 
internal  waves  or  water-mass  boundary  fronts,  as  well  as  absolute  surface  current  measurements 
using  interferometric  SAR.  Fourth,  such  products  could  be  useful  for  near-real-time  location  of 
airmass  boundaries  and  other  phenomena  that  vary  the  boundary  layer  depth.  The  method  may 
also  prove  of  utility  for  analysis  of  sea-surface  clutter  in  the  SPY-1  radar  of  the  AEGIS  system 
aboard  fleet  destroyers  and  cruisers,  allowing  these  ships  to  monitor  the  boundary  layer  depth, 
surface  layer  stability,  and  surface  buoyancy  flux  in  real-time  with  their  tactical  radar.  Finally, 
the  fields  we  provide  might  be  used  as  input  to  radar  propagation  models  (Kerr,  1988;  Babin, 
1995). 
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RELATED  PROJECTS 

This  research  follows  naturally  from  our  work  under  the  HI-RES  ARI  project  that  ended  in  1995. 
The  observations  taken  during  HI-RES,  particularly  the  second  field  deployment,  continue  to 
provide  valuable  data  for  calibration  of  the  nonlinear  convective  structures  model  that  was 
developed  as  part  of  the  HI-RES  project  (Shirer,  et  al.,  1995).  The  footprint  signatures  of  stable 
and  unstable  MABLs  that  provide  the  foundation  for  the  current  work  were  first  characterized  as 
part  of  our  HI-RES  work  (Sikora,  et  al.,  1995;  Sikora,  et  al.,  1997).  In  fact,  it  was  through 
interactions  at  HI-RES  workshops  that  we  developed  the  current  collaborative  project  with  our 
JHUAPL  colleagues,  Robert  Beal  and  Donald  Thompson.  Over  the  long  haul,  this  work  may  have 
a  significant  impact  on  the  ongoing  Lockheed/Penn  State  project  to  develop  techniques  for  using 
the  SPY-1  radar  of  the  AEGIS  destroyers  and  cruisers  to  make  tactically  important 
environmental  measurements  without  disrupting  tactical  utilization  of  the  radar. 
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